Spin-parity dependent tunneling of magnetization in single-molecule magnets 
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Single-molecule magnets facilitate the study of quantum tunneling of magnetization at the meso- 
scopic level. The spin-parity effect is among the fundamental predictions that have yet to be clearly 
observed. It is predicted that quantum tunneling is suppressed at zero transverse field if the total 
spin of the magnetic system is half-integer (Kramers degeneracy) but is allowed in integer spin sys- 
tems. The Landau-Zener method is used to measure the tunnel splitting as a function of transverse 
field. Spin-parity dependent tunneling is established by comparing the transverse field dependence 
of the tunnel splitting of integer and half-integer spin systems. 

PACS numbers: PACS numbers: 75.45. +j, 75.60.Ej 



Single-molecule magnets (SMM) are considered the 
best systems for studying quantum tunneling of mag- 
netization at the mesoscopic level. The first molecule 
shown to be a SMM was Mni2 acetate [[[). It exhibits 
slow magnetization relaxation of its S — 10 ground state 
which is split by axial zero-field splitting. It was the first 
system to show thermally assisted tunneling of magneti- 
zation j| §, [|. Fe 8 and M114 SMMs were the first to 
exhibit ground state tunneling |(| . Tunneling was also 
found in other SMMs (see, for instance, [Q, [| |[). A de- 
tailed study of the influence of environmental degrees of 
freedom on the tunnel process has been started on Fes 
and Mni2 acetate (concerning phonons Jh], [H], nuclear 
spins and dipolar couplings JT^ , |l3| , |l4| ]) which were mo- 
tivated by theories |L5|, [l6|, |17| . The spin-parity effect is 
among the fundamental predictions which have yet to be 
established at the mesoscopic level. It is predicted that 
quantum tunneling is suppressed at zero applied field if 
the total spin of the magnetic system is half-integer but 
is allowed in integer spin systems. Van Hemmen and 
Siito Jll| were the first to suggest the absence of tun- 
neling as a consequence of Kramers degeneracy |]35| . It 
was then shown that tunneling can even be absent with- 
out Kramers degeneracy |n], [2(], ^l). In this case, quan- 
tum phase interference can lead to destructive interfer- 
ence and thus suppress tunneling [2lj| . This effect was 
recently seen in Fe 8 and Mni 2 SMMs H ||). 



There are several reasons why the first attempts || |2J| 
to observe the spin-parity effect were unsuccessful. The 
main reason reflects the influence of environmental de- 
grees of freedom that can induce or suppress tunnel- 
ing: Firstly, hyperfine and dipolar couplings can in- 
duce tunneling via transverse field components; further- 
more, intermolecular superexchange coupling may en- 
hance or suppress tunneling depending on its strength; 
phonons can induce transitions via excited states; and 
finally, faster relaxing species can complicate the inter- 
pretation fbl. 
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FIG. 1: Hysteresis loops of a single crystal of (a) Mn4-(S=9/2) 
and (b) Mn4-(S=8) molecular clusters at different tempera- 
tures and a constant field sweep rate dH z /dt = 0.07 T/s. 



In this letter, we show that these problems can be over- 
come by studying the tunnel splitting as a function of 
transverse field. We selected three SMMs which revealed 
to be sufficiently small to show clearly ground state tun- 
neling. 
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The first is [Mn40 3 (OSiMe 3 )(OAc)3(dbm) 3 ], called 
Mn4-(5 = 9/2), with a half-integer ground state S = 
9/2. The second is [Mn 4 (0 2 CMe) 2 (Hpdm) 6 ][C104]2, 
called Mn 4 -(5 = 8), with an integer ground state S = 
8. The third is the well known Fe 8 SMM with a 5 = 10 
spin ground state ||. The preparation, X-ray struc- 
ture, and detailed physical characterization of both Mn4 
molecules have been presented |)[ p5[ . The complex 
Mii4-(5 = 9/2) crystallizes in a hexagonal space group 
with a crystallographic C3 V symmetry. The complex has 
a trigonal-pyramidal (highly distorted cubane-like) geom- 
etry. This complex is mixed- valent: Mn 3 n Mn Iv . The C3 V 
axis passes through the Mn IV ion and the triply bridg- 
ing siloxide group. DC and AC magnetic susceptibility 
measurements indicate a well isolated S = 9/2 ground 
state J2f|. The complex Mn4-(5 = 8) crystallizes in a 
triclinic lattice. The cation lies on an inversion center 
and consists of a planar Mn4 rhombus that is also mixed- 
valent: MnJpMn!, 1 . DC and AC magnetic susceptibility 
measurements indicate aS = 8 ground state Q. 

All measurements were performed using an array of 
micro-SQUIDs iQ. The high sensitivity of this magne- 
tometer allows us to study single crystals of SMMs of 
sizes of the order of 10 to 500 /im. The field can be 
applied in any direction by separately driving three or- 
thogonal coils. 

Before presenting our measurements we review briefly 
the giant spin model which is the simplest model describ- 
ing the spin system of SMMs. The spin Hamiltonian is 



H = -DSl + iJ trans + smbMo-S • H 



(1) 



S x , S y , and S z are the three components of the spin oper- 
ator; D is the anisotropy constant defining an Ising type 
of anisotropy; H tI& ns, containing S x or S y spin operators, 
gives the transverse anisotropy which is small compared 
to DS Z in SMMs; and the last term describes the Zee- 
man energy associated with an applied field H. This 
Hamiltonian has an energy level spectrum with (25 + 1) 
values which, to a first approximation, can be labeled 
by the quantum numbers m = —5, —(5 — 1), 5 taking 
the z-axis as the quantization axis. The energy spectrum 
can be obtained by using standard diagonalization tech- 
niques. At H = 0, the levels m = ±5 have the lowest 
energy. When a field H z is applied, the levels with m < 
increase in energy, while those with m > decrease. 
Therefore, energy levels of positive and negative quan- 
tum numbers cross at certain values of H z . It turns out 
that the levels cross at fields given by noH z ss nD/gfiB, 
with n = 0, 1, 2, 3, .... 

When the spin Hamiltonian contains transverse terms 
(-Htrans), the level crossings can be "avoided level cross- 
ings". The spin 5 is "in resonance" between two states 
when the local longitudinal field is close to an avoided 
level crossing. The energy gap, the so-called "tunnel 
spitting" A, can be tuned by a transverse field (a field 
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FIG. 2: Field sweeping rate dependence of the tunnel splitting 
A-9/2,9/2 measured by a Landau-Zener method for a M114- 
(S=9/2) samples, at different transverse fields. The Landau- 
Zener method works in the region of high sweeping rates 



where A_ 



9/2,9/2 



is sweeping rate independent. 



applied perpendicular to the S z direction) via the S X H X 
and S y Hy Zeeman terms. 

The effect of these avoided level crossings can be seen 
in hysteresis loop measurements. Figs, la and lb show 
typical hysteresis loop measurements for single crystals of 
Mn 4 -(5 = 9/2) and Mn 4 -(5 = 8) (for Fe 8 data, see f|). 
When the applied field is near an avoided level crossing, 
the magnetization relaxes faster, yielding steps separated 
by plateaus. As the temperature is lowered, there is a 
decrease in the transition rate due to reduced thermally 
assisted tunneling. A similar behavior was observed in 
M1112 acetate clusters [|[ ||[ ||] . As seen for Fes hysteresis 
loops become temperature-independent below 0.4 K in- 
dicating ground state tunneling. The field between two 
resonances allows us to estimate the anisotropy constants 
D. We obtain 0.68 K and 0.43 K for Mn 4 -(5 = 9/2) and 
Mn4-(5 = 8), respectively. 

In order to establish the spin-parity effect, the tunnel 
splitting was measured as a function of transverse field. 
The field dependence of the tunnel splitting is expected 
to be very sensitively dependent on the spin-parity and 
the parity of the avoided level crossing. This approach 
is based on the Landau-Zener model that describes the 
nonadiabatic transition between the two states in a two- 
level system |2^, |27|. The original work by Zener con- 
centrated on the electronic states of a diatomic molecule, 
while Landau considered two atoms that undergo a scat- 
tering process. Their solution of the time-dependent 
Schrodinger equation of a two-level system can be ap- 
plied to many physical systems and it has become an 
important tool for studying tunneling transitions. The 
Landau-Zener model has also been applied to spin tun- 
neling in nanoparticles and clusters 28, j30| , j3l[ j32[ . 



When sweeping the longitudinal field H z at a constant 
rate over an avoided energy level crossing, the tunneling 
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FIG. 3: Tunnel splitting for three SMMs as a function of 
transverse field. The data for Feg where taken from |^2| with 
the transverse field applied along the medium hard axis. 



probability P is given by 



P , — 1 



exp 



ttA 



2 



2hgfiB\m — m'\fiodH z /dt 



(2) 



Here, m and m! are the quantum numbers of the avoided 
level crossing, dH z /dt is the constant field sweep rate, 
g w 2, /iB is the Bohr magneton, and K is Planck's 
constant. 

In order to apply quantitatively the Landau-Zener for- 
mula (Eq. ^]), we first checked the predicted sweep field 
dependence of the tunneling rate. The SMM crystal was 
first placed in a high negative field to saturate the magne- 
tization. Then, the applied field was swept at a constant 
rate over one of the resonance transitions and the frac- 
tion of molecules that reversed their spin was measured. 
The tunnel splitting A was calculated using Eq. ^ and 
is plotted in Fig. 2 as a function of field sweep rate. 
The Landau-Zener method is applicable in the region of 
high sweep rates where A_g/ 2j 9/2 is independent of sweep 
rate . The deviations at lower sweeping rates are mainly 
due to the hole-digging mechanism ]l2] | as observed for 
Fes SMM Jll]] . Such behavior has recently been simu- 
lated || (36f 

Fig. 3 presents the measured tunnel splittings obtained 
by the Landau-Zener method as a function of transverse 
field. Note that for all three SMMs the tunnel splitting 
is finite. However, depending on the spin-parity the sen- 
sitivity to an applied transverse field is completely differ- 
ent. The tunnel splitting increases gradually for an in- 
teger spin, whereas it increases rapidly for a half-integer 
spin. Note also that A is plotted on a logarithmic scale in 
Fig. 3, and that the tunnel probability should depend on 
the second power of A (Eq. ^). Therefore, our measure- 
ments show that the tunneling rate of a half-integer spin 
is strongly transverse field dependent, unlike the case for 



10-7 



<u 
c 
c 
3 



10 




measured 
E = 3.2*10- 2 K 



B. 



'3 4 



1.7*10" 3 K 
3*10" 5 K 



-0.05 0.05 0.1 0.15 
Transverse field (T) 



0.2 



FIG. 4: (a) First attempt to simulated the measured tunnel 
splittings for Mn4-(S I = 9/2). For simplicity, the calculated A 
has been averaged over all possible orientation of the trans- 
verse field, (b) Same as in (a) but taking into account a 
Gaussian distribution of transverse field components with a 
half-width a = 0.035 T. 



an integer spin SMM. 

Figs. 4a and 4b present the first attempt to simu- 
late the measured tunnel splittings for Mn^-^S = 9/2). 
Firstly, one should note that there must be transverse 
terms in the spin Hamiltonian because for -f/trans = 0> 
the measured tunnel splitting should be orders of mag- 
nitude smaller than observed. Secondly, we tried to sim- 
ulate the data with a transverse term (B34/2 * [S Z (S+ + 
S 3 .) + (S^ + S^Sz) with respect to the C 3 symmetry 
of the SMM. However, this term also cannot account for 
the measured tunnel splitting. Finally, we found that ei- 
ther the second order term (E * (S+ + S_)) or a fourth 
order term (B44 * ( S+ + St)) can equally well describe 
the measurements p7| . These results suggest that there 
is a small effect which breaks the C3 symmetry. This 
could be a small strain inside the SMM crystal induced 
by defects. Such defects could result from either a loss of 
solvent molecules or dislocations [ p4[ . 

In the ideal case, one would expect A = for H = 0. 
This is not found in a real SMM because of the influ- 
ence of environmental degrees of freedom that can induce 
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tunneling. In our case, mainly hyperfine and dipolar cou- 
plings induce tunneling via transverse field components. 
Fig. 4b presents a simulation of the measured tunnel 
splitting when taking into account a Gaussian distribu- 
tion of transverse field components with a width a = 
0.035 T. This value is in good agreement with other 
SMMs. 

The simulation of measured tunnel splittings for Mn4- 
(S = 8) is much easier because an integer spin is not 
very sensitive to transverse field components resulting 
from hyperfine and dipolar couplings. We found that a 
second order term with E = 0.057 K can describe well 
the data. 

In conclusion, we have shown that the predicted spin 
parity effect [[l8| [n], |2Cj can indeed be observed by mea- 
suring the tunnel splitting as a function of transverse 
field. An integer spin system is rather insensitive to 
small transverse fields whereas a half-integer spin sys- 
tems is much more sensitive. However, a half-integer spin 
system will also undergo tunneling at zero external field 
as a result of environmental degrees of freedom such as 
hyperfine and dipolar couplings or small intermolecular 
superexchange interaction. 
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[35J The Kramers theorem asserts that no matter how unsym- 
metric the crystal field, an ion possessing an odd number 
of electrons must have a ground state that is at least dou- 
bly degenerate, even in the presence of crystal fields and 
spin-orbit interactions [H. A. Kramers, Proc. Acad. Sci. 
Amsterdam 33, 959 (1930)] The Kramers theorem can 
be found in standard textbooks on quantum mechanics 
L. D. Landau and E. M. Lifschitz, Quantum Mechanics 
(Pergamon, London, 1959). 

[36] A more detailed study show that the tunnel splittings 
obtained by the Landau-Zener method are slightly in- 
fluenced by environnemental effects like hyperfine cou- 
plings |nj . Therefore, one might call it an effective tunnel 
splitting. 

[37] The fit could be improved by using a combination of 
transverse terms but this does not lead to a better in- 
sight concerning the spin parity effect. 



